Priedel oscillations induced surface magnetic anisotropy 
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We present detailed numerical studies of the magnetic anisotropy energy of a magnetic impurity 
near the surface of metallic hosts ( Au and Cu) , that we describe in terms of a realistic tight-binding 
surface Green's function technique. We study the case when spin-orbit coupling originates from the 
d-band of the host material and we also investigate the case of a strong local spin-orbit coupling 
on the impurity itself. The splitting of the impurity's spin-states is calculated to leading order in 
the exchange interaction between the impurity and the host atoms using a diagrammatic Green's 
function technique. The magnetic anisotropy constant is an oscillating function of the separation d 
from the surface: it asymptotically decays as ~ l/d 2 and its oscillation period is determined by the 
extremal vectors of the host's Fermi Surface. Our results clearly show that the host-induced magnetic 
anisotropy energy is by several orders of magnitude smaller than the anisotropy induced by the local 
mechanism, which provides sufficiently large anisotropy values to explain the size dependence of the 
Kondo resistance observed experimentally. 

PACS numbers: 75.20.Hr,75.30.GW 
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I. INTRODUCTION 



It is by now more than fifteen years ago that a sur- 
prising suppression of the Kondo effect in thin films and 
wires of dilute magnetic alloys has been observed 
A few years after the first experiments, Ujsaghy et al. 
proposed that the most likely explanation of the ex- 
perimental observations is a spin-orbit coupling induced 
magnetic anisotropy in the vicinity of the surface of the 
films: In the presence of a surface, spin-orbit (SO) cou- 
pling gives rise to a level splitting of the impurity spin, 
and thus blocks the spin-flip processes responsible for 
the Kondo effecl£&£. Indeed, later experiments seemed 
to be in agreement with this simple scenario and con- 
firmed the predictions that follow from it£ Fitting the 
experimental data for a Au(Fe) film, Ujsaghy et al& es- 
timated the width of the 'dead layer', L c , where the split- 
ting is larger or comparable to the Kondo temperature, 
T K = 0.3-K" ~ 0.03 meV, and obtained L c ~ 180 A. 

To explain the unexpectedly large width of the dead 
layer, Ujsaghy et al. also proposed a model to describe 
surface-anisotropy, which we shall refer to as the host 
spin-orbit coupling (HSO) model. In this model an im- 
purity with a half-filled rf-shell and spin S = 5/2 is im- 
mersed in a host metal, where conduction electrons expe- 
rience SO scattering through hybridizing with low- lying 
valence d-orbitals of the host material^i These calcu- 
lations have been revised recently in Ref. [^. This HSO 
mechanism does not lead to the splitting of the six-fold 
degenerate spin state of the impurity, when placed in a 
bulk host with high (cubic or continuous rotational) sym- 
metry. However, the presence of the surface induces an 



anisotropy term, 



H. 



HSO 



K(d) (nS) 



(1) 



where n is the normal vector of the surface, S is the 
spin-operator, and K(d) denotes the magnetic anisotropy 
constant at a distance d from the surface. The anisotropy 
constant K(d) can be estimated in a simple free electron 
model, by treating the spin-orbit coupling, £, and the 
exchange coupling, J, perturbatively. This calculation 
leads to the asymptotic form^ 



K(d) = A(k F )J 2 £ 



2 sm(2kpd) 
~~d 3 



(2) 



where kp is the Fermi wavenumberi** Unfortunately, the 
constant A(kp) contains some cut-off parameters, which 
make the above formula less predictive for the experi- 
ments. However, ab initio calculations^ indicated that 
this bulk mechanism is too weak to explain the experi- 
mental findings. 

Recently, however, a rather different mechanism has 
been proposed to produce a magnetic anisotropy in the 
vicinity of a surfaced This mechanism, which we shall 
refer to as local spin-orbit coupling (LSO) mechanism, as- 
sumes only a strong local SO coupling on the impurity's 
<i-level. The basic observation leading to this mechanism 
is that, for partially filled d-shclls, spin-states have also a 
large orbital content. Therefore, spin states couple very 
strongly to Friedel oscillations in the vicinity of a surface: 
electrons on the deep d-levels can lower their energy by 
hybridizing with the conduction electrons through vir- 
tual fluctuations. The corresponding anisotropy appears 
already to first order in the exchange coupling J and de- 
cays as ~ l/d 2 . In the specific case of an impurity with 
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a d 1 configuration, the corresponding J 3 / 2 ground state 
multiplet is split by the presence of the surface a*^ 

H^P=K(d)(nJ) 2 , (3) 

where J stands for the total angular momentum oper- 
ators, and K{d) ~ Jsm(Qpd)/d 2 , with Qp being the 
length of an extremal vector of the Fermi Surface (FS). 
As shown in Ref. [H the anisotropy can take the desired 
value of about a few tenths of mcV even beyond 100 A 
from the surface. 

Although the second (local) mechanism is expected to 
be dominant for impurities with partially filled (not half- 
filled) (i-shells, in Ref. [ID only a toy model, namely, a 
single-band metal on a simple cubic lattice, has been con- 
sidered. For a quantitative comparison, however, and to 
decide, which mechanism is responsible for the surface- 
induced anisotropy, more realistic lattice and band struc- 
tures should be used. The aim of the present work is to 
provide such a qualitative and quantitative comparison 
of the two mechanisms described above. For this pur- 
pose, we shall embed the impurity into an fee lattice, 
and employ realistic tight-binding surface Green's func- 
tion methods^ to describe the conduction and valence 
electrons of the host material. This method allows for 
a numerically exact treatment of the surface, and also 
incorporates the SO coupling non-perturbatively. 

To describe the magnetic impurity, we shall integrate 
out virtual charge fluctuations on the d-level of the mag- 
netic impurity, and construct realistic spin models, which 
take into account the specific magnetic and crystal field 
structure of the impurity^ 

We shall then study the surface-induced anisotropy 
within both models, and derive explicit expressions for 
the anisotropy constants in terms of the local density of 
states around the magnetic impurity. Analyzing the be- 
havior of K(d) in the asymptotic regime, we find that the 
oscillations of K(d) are related to the extremal vectors of 
the Fermi surface. For the case of Au and Cu host met- 
als, we perform numerical calculations of the anisotropy 
constants as based on the asymptotic formulas and the 
oscillation periods are directly identified from the numer- 
ically computed Fermi surface. In the case of the local 
SO model, we are also able to confirm numerically the va- 
lidity of the asymptotic expressions. Our results support 
the priority of the local SO mechanism. 

II. SHORT REVIEW OF THE THEORETICAL 
APPROACH 

Before we present our results, let us to some extent out- 
line the theoretical methods we use. As mentioned in the 
introduction, in our approach we describe the host ma- 
terial within a tight binding Green's function formalism. 
The interaction between the magnetic impurity and the 
host, on the other hand, is described in terms of an effec- 
tive interaction, which we construct by combining group 
theoretical methods with many-body techniques. Once 



this effective exchange interaction Hamiltonian at hand, 
we can use relatively standard field theoretical toolsi£ 
to do perturbation theory in the exchange coupling, and 
determine the surface-induced anisotropy. 

A. The Green's function of the host 

In this paper, we shall study surfaces of ordered non- 
magnetic hosts, like the (001) surface of Au or Cu. In 
this case, the Hamiltonian of the host can be written as, 

HZi*> = ( £ a<5aa"W +£(^-s) ActAV ,) V<W + 

+ V™' n '5 aa , , (4) 

where p, p' denote atomic layers normal to the surface, n, 
n' label atomic sites within the layers, while A, A' denote 
the canonical spd orbitals centered at the lattice posi- 
tions and cr, a' are spin indices. In Eq. ([5]), we replaced 
all the parameters by their bulk values, i.e., we neglected 
the dependence of the on-site energies, E\ and the SO 
parameter, £, on the layer depth p. The hopping matrix 
elements, V^^F " are confined to first and next-nearest 
neighbors, and their layer-dependence is also neglected. 
These approximations lead certainly to some errors in 
the calculated electronic structure very close to the sur- 
face, however, they are expected to have no serious con- 
sequences in the asymptotic regime, which is the subject 
of our interest. By the same token, in the vacuum (i.e., 
p < 0) the on-site energies are taken to be infinity. This 
simplifies somewhat our calculations, since only layers 
p > 1, forming thus a semi-infinite system, need to be 
considered in the evaluation of the Green's function. 

The Hamiltonian Eq. ([4]) can be recast into a matrix 
in the spin and orbital labels, 

= • (5) 

Since our system is translational invariant within the lay- 
ers, we can also define the Fourier transform of the Hamil- 

° VP 

tonian matrix, B_ (k), and introduce the 'semi-infinite' 
matrix, 

H(k) = \H PP \k)\ . (6) 

I > p,p'>X 

The resolvent or Green's function matrix is then given as 
usual 

G(z,k)= (z-Ti(k))" 1 , (7) 

with z a complex number (energy). 

To perform the matrix inversion in Eq. ([7]), we used 
the surface Green's function technique^ which is an effi- 
cient and, in principle, exact tool to compute the Green's 
function. Most importantly, the computational time of 
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this method scales linearly with the number of layers, for 
which the Green's function is evaluated. The real-space 
representation of the Green's function can then be ob- 
tained by performing the following Brillouin zone (BZ) 
integral, 

G Pn ' P ' n \z) = J- / d 2 kG PP '(z,k)e^^~^ , (8) 

"BZ J 

where £Ibz is the volume of the 2-dimensional Brillouin 
zone, and the translation vector T ra is related to the po- 
sition of atom n in layer p as R pn = C p + T„, with C p 
a layer-dependent reference position. 

The host-Hamiltonian, Eq.Q, must be slightly modi- 
fied in the presence of a magnetic impurity. In this case, 
the hopping of the conduction electrons to the impurity's 
d-orbitals should be excluded, since these processes in- 
volve charge fluctuations at the magnetic impurity site, 
which will be incorporated in the effective exchange in- 
teraction (see next section) . The simplest way to account 
for this constraint is to shift the on-site <i-state energies 
of the impurity, e\, far below the valence band, and add 
the following term to the Hamiltonian, 

AH \J P x"*' n ' = (4 - ex) 5 Pq 5 plq S n0 <5„<o S x ,\> <W , (9) 

where the impurity is at site n = and in layer q. 

To describe the spin dynamics, we do not need to evalu- 
ate the full Green's function: We need its value only for a 
small cluster of sites, C^ q \ consisting of nearest neighbor 
atoms around the impurity and the impurity itself. For- 
tunately, since AH is also local, this restricted Green's 
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|2> 
|3> 
|4> 

|5) D 2z 2_ x 2_ y 2 



TABLE I: Combinations of s-orbitals centered at the 12 neigh- 
bor sites around an impurity having the symmetry of atomic 
d orbitals. 



To construct the effective interaction between the host 
electrons and the magnetic impurity, one can safely as- 
sume that the deep <i-levels of the magnetic impurity hy- 
bridize only with the s-orbitals of the neighboring host 
atoms. However, by symmetry, the deep <i-levels can 
hybridize only with appropriate linear combinations of 
these s-orbitals, a G {x 2 — y 2 , 2z 2 — x 2 — y 2 , xy, xz, yz}. 
In case of an fee lattice, e.g., we have 12 nearest neigh- 
bor s-orbitals, which can be labeled by s xy , s xy , s xy , s xy , 
. . . , Syz, s yz , the subscripts xy and xy referring to neigh- 
boring sites at the positions a(|,^,0) and a(— 1,0) 



function, 

g} q) { z ) = |g p "- p '"'| c( } , (io) 

can be evaluated as 

g iq) 0) = U- g (9) (z) A^ (9) ) 1 g iq) (z) , (11) 

where I_ is the unit matrix, and the matrix elements of 

g iq) (z) and AH} q) arc defined in Eqs. © and ©, re- 
spectively. Finally, the spectral function matrix on this 
cluster is defined as 

g iq \s) = -^- lim (G (q \e + iS)-G (q) (e-tS)) . 
= 2m 5— >+0 V — — / 

(12) 

As shown in the following subsections, the matrix ele- 
ments of this spectral function matrix arc directly related 
to the magnetic anisotropy. 



B. Host spin-orbit model of the magnetic 
anisotropy 

As in Refs. let us first consider a spin 5 = 5/2 

impurity with a half-filled d-shell. In this case, we can 
neglect the SO interaction on the magnetic ion, and the 
bulk SO interaction is the primary source of the surface- 
induced anisotropy. 



I 

relative to the impurity, respectively, and a denoting the 
cubic lattice constant. However, only 5 out of these 12 
states will have a d-wave character, and hybridize with 
the d-levels of the magnetic impurity. These 5 states are 
listed in Table fl] Using these 5 spin-degenerate states, 
we can perform a Schrieffer- Wolff transformation^ that 
leads to the following Hamiltonian, 

5 

E E J « E i^^^. (13) 

i—x,y,z a — 1 a^a' — ±1 



j 



Dxz — 2 (^xz ~h 5 XZ S X z ^xz) 
Dyz — 2 ( s yz + s yz ~~ s yz ~ s yz) 

D X y — 2 (*^ x y ^xy S X y S x y) 

D x 2_y2 — 2\/2 ^^^ Z ^y z ^y z ^xz Sxz s xz £xz) 

-7= (2S X y -f" 25 X y ~\- 2S X y ~\~ 2S X y — Sy Z Sy Z ^y Z Sy Z S XZ S XZ S XZ SxZ 
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Here s,s' = f denote the z-components of the 

impurity spins, S 1 , and a 1 denote the Pauli matrices. The 
operator c^ aa creates a conduction electron with spin a in 
one of the states \a) listed in Table fl] In the bulk, only 
two of the exchange constants J a are independent, since 



by symmetry we have J ; 



■'■</ 



J, 



J yz and J x 2 



J2z 2 -x 2 -y 2 - bi the following, for the sake of simplicity, we 
shall set all these coupling constants equal, and take J a = 
J. This assumption does not modify our conclusions. 

The anisotropy induced by the surface can be com- 
puted by representing the spin in terms of Abrikosov 
pseudofcrmions, and then doing second order calculation 
in the exchange coupling^ The zero temperature first 
and second order contributions to the static (oj — 0) self- 
energy of the impurity spin can be expressed in terms 
of the local density of states (spectral function) matrix, 



P 



a,a:a a 



^sl= I** deTr{g(e)H JiSS ,} 

J — OO 

= f^Tr^eK} , (14) 



and 



y( 2 ) - 
^ss> — 



de de' 



J — oo J e 



^Tr {Q{e)Hj !sS Q{e')Hj >Ss ,} 

s 

de de' 



e - £ 



Tr{ e (e)a i Q(s')^} , (15) 

with Tr{. . . } denoting the trace in the 10-dimensional 
subspacc of the conduction electrons, and ep the Fermi 
energy. The spectral function, g a .a:a'a', can easily be 
obtained from the real-space spectral function matrix el- 
ements, Eq. fl2]) . 

Exploiting furthermore the tetragonal (C& v ) symmetry 
of an fcc(OOl) surface system and time-reversal invari- 
ance, we find that g a ,a;a'a' has the following structure: 



/ Qih 

-ig 7 a y 



Qih 
iQ6cr y 
-ig 7 a x 



ig e a x ig 7 a y -ig & a y \ 

-igevy igi<7 x ig%<y x 

gih igyiz 

Qi I 2 j 



(16) 



where Vft S M and we dropped the energy argument of 
the spectral functions. The above form of g is fully con- 
firmed by our numerical calculations. Inserting Eq. (I16p 
into Eqs. (14]) and ([15]) yields, E^, = 0, and we find 



,(2) 



CK 



HSO 



{Si 



(17) 



where C is a constant and the anisotropy constant, 
if hso j can be expressed as 



K 



HSO 



^HSO 



-^HSO 



-^HSO 



/v HSO 



K 



HSO i 



(18) 



with 



#hso = -4^ 



oo J Sf 



de I de 



, ft(e)£'i(e') 



e' — e 



(19) 



If the impurity is placed in the bulk, then cubic sym- 
metry further implies that 



Qi{e) = Q2{e) , Qz(e) = gi{e) , g 6 (e) = -g 5 (e) 



<?8 



[e) = V3g 7 {e) , g g (e) = -2g 7 (e) 



(20) 



and we obtain i^HSO = 0. Thus the anisotropy is in- 
deed generated by the surface, which breaks the cubic 
symmetry of the crystal. 



C. The local spin-orbit coupling model of the 
magnetic anisotropy 

As in Ref. IT3 , let us now consider a magnetic impu- 
rity in a d 1 configuration such as a V A+ or Ti 3+ ion. 
In this case, according to Hund's third rule, a strong lo- 
cal spin-orbit coupling will lead to a J = 3/2 multiplct 
that is separated from the J = 5/2 multiplet typically 
by an energy of the order of ~ 1 eV. In a cubic crystal 
field, this J ~ 3/2 ground multiplet remains degener- 
ate (]?8 double representation), implying that no mag- 
netic anisotropy appears if the magnetic impurity is in 
the bulk. Anisotropy will, however, arise, once the im- 
purity is placed to the vicinity of a surface that breaks 
the cubic symmetry. 

To construct the exchange interaction between the con- 
duction electrons and the magnetic impurity, we first 
notice that the impurity's J — 3/2 multiplct can hy- 
bridize only with those linear combinations of neighbor- 
ing s-states which transform according to the same (I^) 
representation. Such a four-dimensional rf-type set can 
be constructed from the states in Tabled] as 



|S-3/2) = D X 2_y2 ||) , 

\ S -l/2) = D 2z 2_ x 2_ y 2 \[) , 

|Sl/2) = D 2z 2_ x 2^ y 2 |t) , 

IS3/2) = -D x 2_ y 2 |t) . 



(21) 
(22) 
(23) 
(24) 



Assuming that the impurity-host interaction is mainly 
dominated by quantum fluctuations to the (non- 
degenerate) d° state, in lowest order of the hybridization, 
a Coqblin-Schrieffer transformation leads to the follow- 
ing effective exchange interaction ; 1 ^ 17 



3 3 

H,j = J s L s m' \-m')(-m\ 



(25) 



where | |m) stand for the four states of the Tg impu- 
rity multiplet, and are creation operators creating an 
electron in the host states ([21]) (Un- 
interestingly, due to the different orbital content of the 
impurity states ||, ±|) and ||,±|), already the first or- 
der contribution to the self-energy gives a non-vanishing 
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anisotropy in the vicinity of a surface^ 

/£ F 
deQmm' (s) 
-oo 



(26) 



The local spectral function of the host is now a 4 x 4 
matrix, p mm '(e), that has a diagonal structure, and is 
related to the spectral functions defined in Eq. p6[) as 
follows, 



Qrava' (^) — Qm (^) &mm' 



(27) 



f?±3/2(e) = Qs{e) , Q±i/2(e) = Qa{e) ■ (28) 

From Eq. (j2TJ)) it is obvious that the J = 3/2 multi- 
plet is degenerate under cubic symmetry (in the bulk), 
while under tetragonal symmetry it is split by an effective 
anisotropy term, Eq. ([3]), with 



LSO 



K 3 

JIT 



LSO -"LSO ) 



and 



K 



LSO 



d£Qi(e) 



(29) 



(30) 



D. Asymptotic form of the anisotropy constants 

The presence of the surface induces Friedel oscillations 
in the local spectral functions^ For large distances d 
from the surface, an asymptotic analysis can be per- 
formed based on the rapid oscillations of the electronic 
wave function, ~ e lfc 2 d _ j n t ne simple case, when the con- 
stant energy surface in the three-dimensional Brillouin 
zone of the bulk system is formed by a single band (like 
the Fermi surface of noble metals), this leads to the fol- 
lowing expressions for the spectral functions appearing 
in Eq. (USD, 

Qi (e, d) ^g° i {s) + - d Y J 9? (e) cos [Q n (e) d + 6? (e)] , 

71 

(31) 

where (e) is the spectral function in the bulk, and the 
Q n (e) 's denote the lengths of extremal vectors of the con- 
stant energy surface, normal to the geometrical surface. 
The factors </" (e) denote the amplitudes of the oscilla- 
tions, and 0™ (e) are their phase. As we shall discuss later, 
in case of an fee (001) geometry there are two different 
extremal vectors. Furthermore, it turns out that each of 
the spectral function matrixelements has a non-negligible 
contribution related only to one of these vectors, there- 
fore, as what follows, we shall label the extremal vectors 
with the index of the matrixelements i. By substituting 
expression (f3Tj) into Eqs. (fl9|) and (|30j) we then obtain 
the asymptotic form of the anisotropy constants. 



1. Host spin-orbit coupling model 

In case of the host spin-orbit coupling model, the con- 
tributions, ATygQ, to the magnetic anisotropy constant 



can be expressed in leading order of l/d as 



4.7 



2 /* oo 



K HSO = d~ Re 1 ~ 



de 



(32) 



CFS 
SF+S 



deg° (e-e) 9i (e) e *N«M«H-**(«0] 



Assuming that p° (e), Qi (e), gi (e) and 6i (e) are slowly 
varying functions of e, whereas e l Qi( e ) d is rapidly os- 
cillating, the inner integrals in Eq. (|3"2"|) give sizable 
contributions only for small values of e, and therefore, 
we can expand Qi{e) around ep, Qi{e) ~ Qi(sp) + 
Qt(sp) (s — ep) and substitute all the other functions by 
their values at Ep. This procedure yields the following 
asymptotic form: 



K 



4 J 2 tt£»° (sp) gi (e f ) cos [Qi (e f ) d + Q % (e F )] 



HSO 



\Q'i {sf) I d 2 

(33) 

For free electrons, Q(ef) — 2kp, and the above re- 
sult resembles that of Ujsaghy et air- however with a 
~ l/d 2 rather than a ~ l/<i 3 decay. This difference is a 
consequence of the assumption made in Ref. that the 
scatterers in the host are distributed homogeneously. 



2. Local spin-orbit coupling model 

In case of the local spin-orbit coupling model the en- 
ergy integral in Eq. (|30p can be easily performed yielding 



K 



J gi (e f ) sin [Qi (e f ) d + { (s F )} 



LSO 



2|Q<(£f)| 



d 2 



(34) 



Interestingly, the asymptotic rf-dependence of the mag- 
netic anisotropy is described by very similar functions 
within both models, only the coefficients and the prcfac- 
tors arc different. 



III. COMPUTATIONAL DETAILS 

For a realistic description of the host's valence and 
conduction bands we used the on-site energies and the 
first and second nearest neighbor hopping parameters as 
given in Ref. [l9| for Au and in Ref. [2(] for Cu, and set 
the cubic lattice constants to their experimental value 
cubic, a C u = 3.615 A and a Au = 4.078 kM- The spin- 
orbit parameter, £, has been determined from the dif- 
ference of the SO-split <i-rcsonance energies, AEd = 
Ej=5/2 ~ Ej—^/2, derived from self-consistent relativis- 
tic first-principles calculations. 22 This splitting is related 
to our SO coupling as 



AE rl 



(35) 



6 



For Au we thus obtained £ = 0.64 eV, while for Cu 
£ = 0.1 cV. In order to reduce the computational ef- 
forts in performing the Brillouin zone integrals, Eq. ((8]), 
we made use of the C^ v point-group symmetry of the 
fcc(OOl) surface and applied an adaptive uniform mesh 
refinement for sampling the fc-points in the irreducible 
(1/8) segment of the Brillouin zone. In general, about 
10 4 fc-points were sufficient to calculate all the spectral 
function matrix elements in (|16p with a relative accuracy 
of 1 %. We performed calculations for the g^s for up 
to 50 monolayers below the surface, corresponding to a 
separation of d ~ 90 A and d ~ 100 A for Cu and Au, 
respectively. 

Performing the double energy integral in Eq. (fT9]) is a 
quite demanding numerical procedure. Therefore, for the 
host spin-orbit model, we first fitted the spectral func- 
tion matrix elements by the function (|3ip . and then used 
the asymptotic form, Eq. (|33[) to compute the magnetic 
anisotropy, AThso- As we shall see later, beyond about 10 
atomic layers (d > 20 A ) the calculated matrix elements 
followed the asymptotic form, and the parameters, <?i(e), 
0i(e) and Qi(e) could be fitted with a high accuracy. 

In case of the local spin-orbit coupling model, we also 
performed a similar procedure to calculate the magnetic 
anisotropy constant in the asymptotic regime, Eq. ([34]) . 
However, in this case, it was also possible to compute 
the anisotropy constant directly from Eq. (|30[) : In this 
case, we could deform the energy integration contour by 
using the analyticity of the Green's function on the com- 
plex plane, and as few as 12 energy points along a semi- 
circular contour in the upper complex half-plane were suf- 
ficient for a very accurate evaluation of the corresponding 
integral. 



IV. RESULTS 



A. Electronic structure of the bulk host 

We first performed calculations of the density of states 
(DOS) of bulk Cu and Au. As shown in Fig. [TJ the 
dispersion of the 3c?-band of Cu is about 4 eV, while 
the 5d-band of Au is much broader (~ 7 cV). Reassur- 
ingly, the positions and the heights of the characteristic 
peaks of the DOS compare well with those obtained from 
self-consistent first principles calculations^^ Clearly, in 
copper, the small SO coupling, £ = 0.1 eV, causes merely 
a slight modification in the DOS in the vicinity of the d- 
like on-site energy (~ 5.07 eV). In the case of Au the 
SO coupling is much stronger, £ = 0.64 eV, and is large 
enough to influence the whole d-band: It gives rise to 
strong splittings of the dispersion peaks and it also in- 
creases slightly the bandwidth. As indicated by the ver- 
tical lines in Fig.[TJ the Fermi energies, Ep u = 8.3 eV and 
£p u = 7.4 cV, lie well above the (i-band for both metals. 

As we learned from the asymptotic analysis presented 
in Sec. II. D, extremal vectors of the Fermi surface play a 
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FIG. 1: (Color online) Calculated valence band densities of 
states for Cu and Au bulk without SO interaction (dots) and 
with SO interaction (solid line) . For the latter case the Fermi 
energies, Ep u = 8.3 eV and Ep u = 7.4 eV, are indicated by 
vertical lines. 
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FIG. 2: (Color online) Calculated plane cuts perpendic- 
ular to the (1 -1 0) direction of the FS of Cu and Au. 
The arrows denote the extremal vectors of lengths, Q^ n = 
0.505 A" 1 , Q^ax = 1-208 A" 1 and Q*$ n = 0.298 A" 1 , 



'max 



= 1.228 A -1 . 



crucial role in determining the magnetic anisotropy con- 
stants. Therefore, we next investigated the plane cuts of 
the Fermi surface perpendicular to the (1 -1 0) direction. 
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One can easily read off the length of the (001) extremal 
vectors from the cuts depicted in Fig. [5] The absolute 
minimum of the width of the Fermi surface, Q m i n , can 
be found at k = 0, while the maximum width of the cor- 
responding cut, Qmax, is related to saddle-points of the 
Fermi surface. In the case of a Cu host the values ob- 
tained from our tight-binding analysis, Q^ n = 0.505 A -1 
and <2^ax = I- 208 A -1 correspond to periods of 12.44 A 
and 5.20 A (6.88 and 2.88 monolayers (ML)) of the oscil- 
lations, and agree fairly well with the periods, 6.08 ML 
and 2.60 ML, calculated by Lathiotakis et ali^. Similar 
satisfactory agreement can be found in the case of a Au 
host between the periods found from our present calcu- 
lations, 10.34 ML and 2.51 ML, and those calculated by 
Bruno and Chappert, 8.6 ML and 2.6 ML.^ It should be 
noted, however, that the shape of the FS depends very 
sensitively on the position of the Fermi energy the precise 
determination of which is a quite subtle task, since for 
noble metals like Cu and Au the Fermi energy lies in the 
very flat 4sp band (see also Fig. []}. 



The magnetic anisotropy constants within the 
host spin-orbit coupling model 



B. 



We calculated the spectral function matrices, Eq. ([161) , 
at the Fermi energy of Cu and Au using the methods 
described in Sections II. B. and C, for up to 50 ML be- 
low the surface. As a convincing check of our numeri- 
cal procedure wc verified that the structure of the cal- 
culated matrices agrees with that derived analytically 
from symmetry principles. In the case of a Au host, in 
Fig. [3] we plotted the calculated off-diagonal matrix ele- 
ments, £>5(£f), • ■ ■ , <?9(£f), as a function of the distance 
d from the surface. As expected, large oscillations can be 
observed for all the spectral functions near the surface 
(d < 20 A). These oscillations, however, survive for large 
distances only for £?g , while they are strongly damped 
in all the other cases. The limiting values of Qi corre- 
spond to the bulk case and, as we checked, satisfy the 
conditions, Eq. (j2"0)) with less than 1% relative numerical 
accuracy. 

In Fig. |4] we display the spectral function g§(e F ) on 
an enlarged scale, together with a fitting function of the 
form, Eq. (|31[) . Quite surprisingly, the asymptotic func- 
tion applies even in the range of d > 20 A and, therefore, 
there is no need to perform a 'preasymptotic' analysis as 
suggested in Ref. [t| The fitted parameters of Eq. ([3T]) arc 
as follows: g%{e F ) = -3.99 ± 0.01 • 10~ 4 eV _1 , g e (e F ) = 
-1.484 ± 0.008 • 10" 3 AcV" 1 , Q 6 (e F ) = 1.2228 ± 0.0001 
A" 1 , and 9 6 (e F ) = 1.324 ± 0.006 rad. It is particularly 
noteworthy that the fitted wavenumber agrees with an 
accuracy of 0.5 % with the length of the extremal vec- 
tor, Qmax, computed from the Au Fermi Surface. We 
could fit all other off-diagonal spectral function compo- 
nents entering the expression of -f^HSO with a similar fit 
with exactly the same wavenumber. However, the am- 
plitude of these other components was by at least two 




100 
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FIG. 3: (Color online) Calculated off-diagonal spectral func- 
tion matrix elements (see Eq. (|16[) h at the Fermi energy as 
a function of the distance, d, from the (001) surface of a Au 
host. 
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FIG. 4: (Color online) Asymptotic fit to the function pip 
(solid line) of the calculated values of the Qs(ef) spectral func- 
tion (triangles) as a function of the distance from a Au(001) 
surface. The dashed line denotes the bulk value of q§(ef)- 



orders of magnitude smaller than g§(e F ). 

Our calculations thus indicate that the long- 
wavelength oscillation corresponding to Q m i n of the FS 
either enters with a negligibly small amplitude or doesn't 
enter at all in the asymptotic form of the off-diagonal 
spectral function matrixclcmcnts. This can easily be 
understood by noticing that the asymptotic contribu- 
tions to the real-space spectral function matrixclcmcnts, 
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d (A) 

FIG. 5: (Color online) Upper panel: The K% so contribution 
to the magnetic anisotropy constant within the host spin- 
orbit coupling model for a Au host as calculated from the 
asymptotic expression, Eq. (|33|) . Lower panel: The K% so 
contribution to the magnetic anisotropy constant in the case 
of a Cu host. In both cases an exchange interaction parameter 
J = 1 eV was used. 



above £f and obtained Q'(ef) = 0.235 (A cV)" 1 . Thus, 
in case of a Au host we get the following asymptotic 
function for Kfj SO (d) (displayed in the upper panel of 
Fig. ED 

K Hso( d ) = cos [1-2228 ■ d + 1.324] ficV , (37) 

where d is measured in A. Notice the surprisingly small 
magnitude of Kfj SO : even at a distance of about d = 20 
A the amplitude of the above oscillating function is about 
0.079 ficV. 

Wc performed similar calculations for a Cu host. In 
Cu, the spectral functions show asymptotic oscillations 
with Q(sf) = 1-205 A" 1 that agrees within 0.3 % with 
the length of the extremal vector, Q max , of the Cu FS. 
In Cu, the Kfjso contribution shown in the lower panel 
of Fig. [5] dominates the magnetic anisotropy. This is in 
the range of 0.01 neV = 10~ n cV, i.e., it is at least by 
three orders of magnitude smaller than the one found 
in case a Au host. This decrease is mostly due to the 
smaller SO interaction in Cu than in Au. As we checked 
by varying £ for Au, the spectral functions in Eq. (fl~6|) 
scale linearly with £, therefore, by Eq. (fP9|) the magnetic 
anisotropy constant scales as ~ £ 2 . This result clearly 
justifies the approach of Ujsaghy et a/., who treated the 
SO interaction perturbatively.^^. 



^+ P K( ?+P ')«' (£) {p p) = Q> d = q a j rdatcd tQ Q m . n 

are of the following form, 

Pip' i \ 

(q+p)n,(q+p')n' /\ ^ (0)p n,p' n' , s 9 s <j,s<j'\ £ ! 
"sa.sa' V°/ — tser.ser' \r) d 

cos[Q min (e)d + 9{e)} , (36) 

where gf^ p s ^! p ™ (e) refer to the corresponding bulk ma- 
trixelements. Eq. (|36[) implies that the oscillating part 
does not depend on the in-plane positions, n and n', 
which is the consequence that the minimal extremal vec- 
tor is at the k = point of the 2D Brillouin zone. As 
explained in Sec. II. B., the matrixelements in Eq. (|16p 
are linear combinations of the above real-space matrix 
elements according to the states in Table I. Since the 
states |a) (a = 1, ... ,4) are constructed as antisymmet- 
ric combinations of neighboring s-orbitals in the same 
plane, q + p, or as a sum of such antisymmetric combi- 
nations, in their matrixelements the asymptotic oscilla- 
tory part corresponding to Q m i n necessarily cancels. As a 
consequence, only the spectral function g^ = (5\g\5) has 
asymptotic oscillations with wavenumber Q m [ n , which, 
however, does not give a contribution in the host SO 
model. 

We calculated the magnetic anisotropy constant using 
the asymptotic fits of the spectral functions and Eq. (f33|) . 
We numerically determined the energy derivative of the 
magnitude of the extremal vector, Q'(sf), by fitting the 
spectral functions at two energy values close below and 



C. The magnetic anisotropy constants within the 
local spin-orbit coupling model 

As pointed out in Ref. [Tl, a mechanism based on a 
strong local SO interaction of the impurity (local SO 
model) can give rise to a level splitting that is orders of 
magnitude larger than the host-induced anisotropy. To 
demonstrate this idea, in Ref. [l2| we studied the simple 
but unrealistic case of a single-band metal on a simple 
cubic lattice. Here wc extend the calculations of Ref. LL2| 
and perform calculations for realistic host metals (Cu and 
Au). 

According to the theory presented in Sec. II. C, we 
need to compute the diagonal spectral function matrix- 
elements, g 3 = (4|g|4) and £4 = (5|p|5), see Table I and 
Eq. (fTB"|) . Our calculations clearly showed that the d de- 
pendent Fricdel oscillations of g% are several order smaller 
in magnitude than those of £14. This can be understood 
by noticing that, due to the different spatial character of 
these two states {D x 2_ y i and D 2z 2_ x 2_ y 2), g 3 comprises 
an average of spectral weights in layers q — 1 and q + 1 , 
while g>4 takes the difference of spectral weights in layer 
q with respect those in layers q — 1 and q + 1, q denoting 
the layer of the impurity's position. Recalling that for 
a cubic bulk £3 = ^4, see Eq. (|20p . in the asymptotic 
region i^Lso becomes proportional with the integral of 
the function, Ag4(e, d) = g^e.d) — Q%{e). This function 
is displayed in Fig. [S]for both the Au and the Cu hosts. 
Remarkably, the amplitude of the Fricdel oscillations is 
about one order of magnitude larger than those of the 



9 



0.20- 



0.00 



-0.20 



1.00 
0.50 
0.00 
-0.50 
-1.00 







■ i 


, 1 1 1 1 1 1 

Cu 

A/vw 


\ ' i 


1 1 1 1 1 1 1 1 1 1 1 1 1 



FIG. 6 

Eq, 



20 30 40 50 60 70 80 90 100 
d (A) 

(Color online) Calculated values of Aq4(ef) = 
Qa{sf) (squares) with corresponding asymptotic fits, 
, (solid line) as a function of the distance from the 




(001) surface of Au and Cu. 
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FIG. 7: (Color online) Magnetic anisotropy constants within 
the local spin-orbit coupling model calculated by using the 
asymptotic formula, Eq. (|34[) . as a function of the distance 
d from the (001) surface of a Au (dashes) and a Cu host 
(solid line). In case of Au the squares stand for the magnetic 
anisotropy constants calculated directly from Eq. (|30p . 



SUMMARY AND CONCLUSIONS 



off-diagonal spectral functions (compare with Fig. [3] for 
the case of Au). Note that the off-diagonal matrixele- 
ments appear in first order of the spin-orbit coupling. 
The oscillations have larger periods as compared to the 
off-diagonal spectral functions: a fit to the asymptotic 
function, Eq. (j3Tj) . shown also in Fig. [6] gave the val- 
ues Q Au = 0.292 A^ 1 and Q Cu = 0.505 A" 1 , in very 
good agreement with the length of the small extremal 
vector, Qmin, of the corresponding Fermi surfaces. Inter- 
estingly, the amplitude of the oscillations is more than 
three times larger for Cu than for Au: from the fits we 
obtained g A (e F ) = 1.16 • 10~ 2 A cV" 1 and 3.53 • 10~ 2 A 
cV _1 for the case of Au and Cu, respectively. 

Fig. [7] shows the magnetic anisotropy constants ob- 
tained using Eq. (fM|) with the parameters extracted from 
the fits of Aq4(ef, d). The parameter Q'(cf) was com- 
puted as for the off-diagonal spectral functions, and took 
a value of 0.245 (AeV)" 1 for Au and 0.238 (AeV)" 1 for 
Cu. Choosing again J = 1 eV, we obtained for the am- 
plitudes of the oscillations of K, A(d) = 0.0237/d 2 eV 
and A(d) = 0.0742/g? 2 eV (d measured in A ) in Au and 
Cu, respectively. In particular, for Cu this gives an am- 
plitude of 0.03 meV at d = 50 A or 0.007 meV at d = 100 
A, which is in the range of Tk for typical dilute magnetic 
alloys, such as Cu(Mn), Cu(Cr). 

In Fig. we also compare the magnetic anisotropy 
constants obtained from the asymptotic analysis with the 
values we get by performing the contour integration in 
Eq. (f3TJ)) . Apparently, already for d > 35 A, these val- 
ues lie almost perfectly on the asymptotic curve. This 
nice agreement proves the validity of the asymptotic for- 
mula, Eq. (|34[) . as well as the accuracy of our numerical 
procedure to compute the magnetic anisotropy constant. 



In this paper, we performed a theoretical study of 
two mechanisms for surface-induced magnetic anisotropy 
of a magnetic impurity: a local spin-orbit mechanism 
(LSO),— and a host spin-orbit mechanism (HSO).— Both 
mechanisms appear as a result of Friedel-like oscillations 
in the local spectral functions, induced by the surface. In 
the local SO mechanism, the rather large diagonal, i.e., 
charge oscillations couple through the local spin-orbit 
coupling on the d- or f-levcl of the magnetic impurity 
to the impurity spin and lead to a surface-induced split- 
ting of the spin states. The host SO mechanism, on the 
other hand, relies on oscillations in the off-diagonal ele- 
ments of the local spectral functions, i.e., oscillations in 
the "spin sector" that couple directly to the spin through 
an exchange interaction. These oscillations are induced 
by the SO coupling in the host metal and, thus, they are 
much weaker than the Friedel oscillations in the "charge 
sector". Based upon this simple picture, one therefore 
expects that the first mechanism is dominant for impuri- 
ties with a partially filled d- or f-shell, while the host SO 
mechanism may become important for half-filled shells, 
in which case the local SO mechanism cannot be at work. 

In this paper we attempted to compare these two 
mechanisms quantitatively. For the description of the 
host's valence and conduction electrons we used a tight- 
binding Green's function technique, which allows for a 
perfect treatment of the semi-infinite surface geometry, 
and makes also possible a non-perturbative treatment of 
the host SO interaction. We then used a field theoretical 
approach to compute the self-energy of the spin up to first 
(local SO) and second orders in the exchange coupling, J, 
(host SO model) , and derived explicit expressions for the 
anisotropy constants, K, as a function of the separation 
d between the impurity and the surface. 
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These expressions have been analyzed using an asymp- 
totic analysis which resulted in a very similar oscillatory 
dependence of K on d in both models: the periods of the 
oscillations could be identified as the magnitudes of the 
extremal vectors of the Fermi Surface of the bulk host 
and their amplitudes decayed in both models as l/d 2 . 
Here we must remark that in our calculations in Ref. [12 
we predicted a 1 /d 3 decay of the oscillations of K within 
the host spin-orbit mechanism. This must be contrasted 
to the results of the present work, where we find rather 
a l/d 2 scaling of the host-induced anisotropy. This ap- 
parent controversy is due to a small difference in the cal- 
culations: Unlike the present work, in Ref. [IH we ne- 
glected the potential scattering at the impurity site, i.e., 
we used the local spectral functions of a perfect semi- 
infinite host. In this case, however, one can show that 
certain off-diagonal elements of the local spectral func- 
tion matrix must vanish due to two-dimensional trans- 
lational symmetry. These off-diagonal matrix elements 
are non-zero, once translational invariance is broken by 
potential scattering at the impurity site, and they give 
rise to a l/d 2 decay of the anisotropy as shown in Sec. 
II.D. 

Using realistic tight-binding parameters, we calculated 
the amplitudes of the magnetic anisotropy oscillations for 
the cases of Au and Cu metal hosts. As expected from 
the very different SO interactions in these metals, within 
the host SO model, the magnetic anisotropy constant for 
Au turned out to be about three orders larger in magni- 
tude than for Cu. Nevertheless, even for a Au host and 
close to the surface, the magnetic anisotropy constants 
remained below the range of 0.1 fj,cV. Though a direct 



comparison with the result of Ref. is quite questionable 
mainly due to the different geometrical distribution of 
the host atoms and to the different approximations used, 
the above value is close to the lower limit of the esti- 
mated range of K given in Ref. 0. We therefore conclude 
that most probably the host SO mechanism of Ref. is 
too weak to explain the size-dependence of the Kondo 
resistance. 

The local SO mechanism proposed in Ref. [IH, on the 
other hand, gives a magnetic anisotropy constants for Cu 
in the range of 0.03-0.01 meV for even at distances 50- 
100 A away from the surface. Although they are in the 
same range, the magnetic anisotropy constants for Au 
were about three times smaller than the ones we got for 
Cu. Our numerical studies imply that the primary mech- 
anism to produce a magnetic anisotropy in the vicinity of 
a surface is provided by the local SO coupling, where the 
local Hund's rule coupling conspires with Friedel oscilla- 
tions to produce a large anisotropy effect. This mecha- 
nism seems to be large enough to explain the suppression 
of the Kondo resistance anomaly observed in thin films 
and it also supposed to be the dominant source of (ran- 
dom) magnetic anisotropy in metallic mcsoscopic struc- 
tures such as metallic nano-grains, nano-wires, or point 
contacts. 
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nanced by the Hungarian National Scientific Research 
Foundation (OTKA T068312 and NF061726) and by 
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and the Hungarian Science and Technology Foundation 
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